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ABSTRACT 

Close-in cxtrasolar gas giants - the hot Jupiters - display departures in radius above the zero- 
tempcrature solution, the radius excess, that are anomalously high. The radius excess of hot Jupiters 
follows a relatively close relation with thermal tidal tidal torques and holds for ^4 — 5 orders of 
magnitude in a characteristic thermal tidal power in such a way that is consistent with basic theoretical 
expectations. The relation suggests that thermal tidal torques determine the global thermodynamic 
and spin state of the hot Jupiters. On empirical grounds, it is shown that theories of hot Jupiter 
inflation that invoke a constant fraction of the stellar flux to be deposited at great depth are, essentially, 
falsified. 
Subject headings: planets: hot Jupiters 
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1. INTRODUCTION 

The radii of hot Jupiters are anomolously larger than 
their zero temperature radii, i?o (Bodenheimer et al. 
2003; Arras & Bildsten 2007). The powerful stellar ra- 
diation field serves as an ample source of energy which 
may inflate hot Jupiters to their current radii (Guillot & 
Showman 2002). The major theoretical challenge is in 
finding a mechanism that can transfer a small amount 
(< 1%) of the stellar insolation to great depth. Current 
ideas include thermal tides (Arras & Socrates 2009a, b; 
2010), delayed gravitational contraction (Ibigui et al. 
2010) and ohmic dissipation due to magnetic activity 
on the heavily irradiated surface (Batygin & Stevenson 
2010). 

It seems that there is an relation between observed 
radius excess A defined as 



A^ Rp 
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Rp 

where R p is the measured planet radius, and a charac- 
teristic thermal tidal power L t hT defined as 
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where cr SB is the Stefan-Bolt zmann constant, c p is the 
specific heat, n is the mean motion and T eq is its equi- 
librium surface temperature. 

The radius excess A is determined by a balance of heat- 
ing and cooling. In the thermal tide scenario, the ther- 
mal tide torque is balanced by the gravitational torque. 
While the thermal tide torque is non-dissipative, the 
gravitational tide torque necessarily is the result of dis- 
sipation. Torque balance in steady state then implies 
continuous heating as a result of gravitational tidal dis- 
sipation. It follows that the planet radius is determined 
by a balance of this steady state heating with cooling at 
the photosphere. The ultimate source of energy is the 
stellar radiation field, which performs work by moving 
material across the tidal potential, near the surface of 
the irradiated planet (Arras & Socrates 2009a; 2010). 

The plan of this paper is as follows: In $2] we motivate 
the need for thermal tidal torques in hot Jupiters, com- 
pute their characteristic values for both equilibrium and 



dynamical thermal tides as well as the the characteristic 
thermal tidal power, all in terms of physical quantities 
that can be observationally inferred. The thermal tide 
- radius excess relation is presented in Sj3j A brief dis- 
cussion is contained in 2J A brief summary is given in 

m 

2. THEORY OF THERMAL TIDES IN HOT JUPITERS 

The radiation field of a star forces the atmospheres of 
planets in orbit. The response of a planetary atmosphere 
to this forcing is referred to as the "thermal tide." Planets 
that spin asynchronously experience a sharp change in 
thermal forcing during twilight hours and consequently, 
a broad spectrum of tidal harmonics are excited. For 
circular orbits, the semi-diurnal (I = m = 2) harmonic 
of the forcing leads to a thermal tidal bulge that can 
couple to the leading order term of the gravitational tidal 
potential. The atmosphere's finite thermal inertia leads 
to a lag in the response with respect to the forcing, which 
allows the tidal potential to exert a net torque on the 
thermal tide during the heating cycle. Since the thermal 
tidal bulge peaks in the morning and after sunset, but 
before midnight, the acceleration increases the planet's 
rotational energy. 

The picture described above has a long history. On 
the basis of the measured magnitude and phase of Earth's 
thermal tide, Thompson (Lord Kelvin) computed Earth's 
thermal tidal torque and found it to be approximately 
equal to one tenth the dissipative lunar tidal torque 
(Thompson 1882; see also Munk & Macdonald 1960). 
Due to its heavy atmosphere and correspondingly large 
thermal tides, Gold & Soter (1969) hypothesized that 
the current spin state of Venus results from a balance of 
thermal and gravitational tidal torques so that its cur- 
rent state of spin can be maintained over the age of the 
solar system. 

In what follows, we motivate the importance of thermal 
tides and thermal tidal torques in hot Jupiters, within 
the context of understanding their inflated radii. 

2.1. tidal power as a source of dissipation 

The characteristic value for the maximum instanta- 
neous tidal power for a uniform spherical body of radius 
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Fig. 1. — The A — LthT relation for a sample of 63 hot Jupitcrs 
with Msini > 0.5 Mj. An approximation for the zero temperature 
radius Rq is obtained by using the cooling models of Baraffe et al. 
(2010) for core-less planets with Z = 0.02 in the absence of irradi- 
ation at an age of 7 Gyrs. Large red points are theory, computed 
using the Gold-Soter approximation for a Jupiter-mass, core-less 
Z = 0.02 planet whose optical absorption opacity k* = 0.003 (Ar- 
ras & Socrates 2009a). The tidal quality factor of the planet is 
chosen such that Q = 10 5 for the Jupiter-Io interaction and the 
constant lag time model of Hut (1981) is employed in order to 
compute the gravitational tidal torque and equilibrium dissipation 
rate. 



2.2. the role of thermal tides 

Stellar irradiation creates a tidal bulge that leads max- 
imum forcing at noon, where the phase shift results from 
thermal inertia. The tidal gravitational field of the star 
torques the thermal tide, accelerating the planet away 
from a state of synchronous spin. As the rate of spin 
increases, the amount of energy absorbed per cycle de- 
creases and consequently, so does the magnitude of the 
thermal tidal quadrupolc. The thermal tide torque even- 
tually comes into balance with the usual dissipative grav- 
itational tidal spin-down torque. In steady-state, the 
planet's spin is asynchronous and therefore, gravitational 
tidal power is continuously dissipated. The ultimate 
source of energy is the starlight of the primary, which 
performs work by moving material across the tidal po- 
tential (Arras & Socrates 2009a; 2010). 

A fundamental underlying assumption of this picture 
is that dissipation of the gravitational tide takes place at 
great depth, where the pressure scale height is compara- 
ble to the radius of the planet. 

Also note that, in the absence of thermal tidal toques, 
gravitational tidal torques are likely to synchronize the 
planet spin extremely fast if the relative strength of tidal 
dissipation in hot Jupiters is comparable to what is com- 
monly inferred from the Jupiter-Io interaction (Goldrcich 
& Soter 1966; Socrates et al. 2012). 



R in a Keplerian orbit about a central source of mass M 
at a heliocentric distance D is given by 
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G 



where the characteristic tidal speed c T is given by 
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and G is Newton's constant 1 . The expression above is 
obtained by taking the peak energy of the equilibrium 
tidal response to dissipate with a timescale given by the 
inverse of y/GM/D s . 

Most hot Jupitcrs are in nearly circular orbits. In that 
case, n = y^GM/a 3 is the mean motion and only the syn- 
chronization tide can lead to tidal dissipation and there- 
fore, the forcing frequency is given by the difference of 
the spin frequency and mean motion. If the planet is 
rapidly spinning such that the characteristic tidal forc- 
ing frequency in eq. [I] is taken to be the mean motion n, 
then eq. Q] may be written as 

-5 



imaxT ~ 10 P; 



r?5 



(3) 



where P3 is the orbital period P or b in units of 3 days and 
Rio is the planet's radius, R p , normalized to 10 10 cm. 
The power required to maintain the radii of heavily in- 
flated hot Jupiters is of order ps 10 28 erg/s and therefore, 
only a small fraction of the total available tidal power is 
required to produce inflated hot Jupiters in steady-state. 

1 Note the similarity of eq. [TJwith maximum luminosity of any 
object in the Universe 



2.3. spin equilibrum and steady tidal power: Gold-Soter 
approximation 

Balance between the gravitational tidal torque and the 
thermal tidal torque determines the planet's spin rate 
and therefore, determines the tidal forcing frequency ui. 
For a circular or nearly circular orbit, torque balance is 
obtained by equating the respective quadrupoles induced 
by thermal and gravitational forcing. The portion of the 
gravitationally excited quadrupole responsible for dissi- 
pation and secular evolution is given by 



n 2 Rl 

Qgvav ~ ~^ 
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(4) 



where the above expression for the tidal quality factor 
Q reflects the frequency dependence of the constant lag 
time model of Hut (1981). Here Qj and wj is the tidal 
quality factor(i=s 10 5 ; Goldreich & Soter 1966) and forc- 
ing frequency (~ 10 hours) of the Jupiter-Io interac- 
tion. Note that we are, equivalently, setting the lag time 
t ss 0.1 s, consistent with the resonant configuration of 
the Galilean satellites (Socrates et al. 2012; Leconte et 
al. 2010; see also Socrates & Katz 2012). 

The Gold-Soter quadrupole due to thermal forcing is 
aproximatcly given by 
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where AM and ith are the mass and thermal relaxation 
time, respectively, of the absorbing layer and ith- In 
terms of physical quantities, the the thermal time may 
be written as 
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By equating eq. [4] with eq. OH the equilibrium forcing 
frequency becomes 
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where T2000J -P4, -RiOi Qs are the equilibrium temperature 
in units 2000 K, the orbital period normalized to four 
days, dimensionless planet radius normalized to 10 10 cm 
and the tidal quality factor of the Jupiter- Io interaction, 
normalized to 10 5 . 

In the equilibrium - Gold-Soter - approximation, the 
tidal dissipation rate resulting from the balance of ther- 
mal and gravitational tidal torques, which we refer to as 
the thermal tidal power L?^ T is given by 
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(8) 



2.4. spin equilibrium and steady tidal power: dynamical 
thermal tides 

Consider the computed thermal tide quadrupolc from 
Arras and Socrates (2010; see their figure 5) for the case 
with outgoing radiation boundary condition. We imag- 
ine an evolutionary scenario in which the planet was born 
spinning rapidly and was then placed in its current orbit 
via some migration mechanism. As it spins down, the 
forcing frequency approaches the g-mode "bump" in the 
thermal tidal quadrupolc, which is responsible for arrest- 
ing any further spin evolution. The first spin up feature 
of the g-mode bump is located near the cutoff frequency 

LO 



TT 



(9) 



where H is scale height of the isothermal layer that lies 
above the convection zone and N is the Brunt- Vaisalla 
frequency, the characteristic rate at which the buoyancy 
force restores motion. 

It is worth noting that the cutoff frequency lu$ is close 
to the mean motion n of hot Jupiters. If wo <C n, the 
quadrupole moments associated with low radial order 
gravity waves would be small and therefore, so would 
the quadrupolc of the thermally-forced response. From 
this perspective, it is somewhat of a co-incidence that 
dynamical thermal tides in the hot Jupiters can lead to 
significant spin-up tidal torques. 

The spin-up quadrupole due to the first feature in the 
3- mode bump may be expressed as 



1th = 9th (u = Uo)f(u,uo) 



(10) 



where / can be approximated by a sharply-peaked Gaus- 
sian of the form 



f(w,iJo) ~ (3e 



_( w _ W0 ) 2 /r 2 



(11) 



where /3 w 0.1 quantifies the amount by which gravity 
waves near the cut-off frequency do not perfectly overlap 



with both thermal and gravitational forcing. The width 
of the Gaussian is narrow so that T <C ujq. 

Spin equilibrium is determined by setting q^ h = q sr av, 
which gives 
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The right hand side must be less than unity in order for 
torque balance to be satisfied. For the physical parame- 
ters that describe the hot Jupiter population, torque bal- 
ance is only possible for Q, > 10 5 . If there is some level 
of sub-surface reflection the excited gravity waves are 
quantized and their response is therefore, characterized 
by a Lorentzian profile near resonance, implying that the 
induced quadrupolar response can be much larger. Con- 
sequently, values of Q 3 ps 10 4 or r ~ 1 s are then possi- 
ble for a steady-state spin equilibrium (Arras & Socrates 
2010). 

Due to the sharpness of the Gaussian, the forcing fre- 
quency in spin equilibrium, given by eq. Q21 is pinned to 
a value that is close to the cutoff frequency 

u ~ Wo- (13) 

By following the steps in eq. |S1 a similar expression for 
thermal tidal power, LJF^, is derived due to the presence 
of dynamical tides 

L*& « 1.2 x 10 



28, 2 p-4 
W 2 -T 4 



R%Q^ erg/s 



(14) 



where w 2 and Q5 are the semi-diurnal forcing frequency 
normalized to (2ir/2 days) and the tidal quality factor of 
the Jupiter- Io interaction, normalized to 10 5 . 

The thermal tidal dissipation rate, eq. [14] resulting 
from the dynamical excitation of gravity waves, is similar 
in form and magnitude to that resulting from the Gold- 
Soter approximation, given by eq. [8] From here on, 
when comparing with data, we employ the Gold-Soter 
approximation since (i) it is simple (ii) in terms of spin 
equilibrium and persistent tidal luminosity, it roughly 
agrees with the more realistic dynamical tide results of 
Arras & Socrates (2010) (iii) the theory of dynamical 
thermal tides is currently in a state of development. 

3. THERMAL TIDE-RADIUS EXCESS RELATION 
The fractional radius anomaly A, given by 



A = 1 - R /R 



v 



(15) 

where Rq is a fiducial radius that is close the zero tem- 
perature solution depends upon planet mass and com- 
position. The value for Rq is obtained by using the on- 
line table provided by Baraffe et al. (2008) by choosing 
choosing Z = 0.02 with no irradiation and their largest 
computed age of 7 Gyr. 

The primary reason for utilizing the models of Baraffe 
et al. for the values of Ro, rather than the actual zero 
temperature or isothermal radius is that they are publicly 
available. Their solutions serve as a fiducial value of 
low planetary entropy that depends on planet mass. It 
follows that any observed departure in planet radius R p 
above Ro necessitates a source of thermal energy from 
within the planet. 

The sample of hot Jupiters is taken from the Exoplanct 
Data Explorer, located on the world wide web. Only hot 
Jupiters with measured mass, radius, effective temper- 
ature and orbit that have measurements i.e., a sizable 



fraction of the total sample, are considered. Further- 
more, the sample is limited to hot Jupiters with masses 
M > Q.5Mj. All of the planets considered have both 
transit and radial velocity measurements, for a total of 
63 objects. 

In figure [1] we plot radius anomaly A against the dis- 
sipation rate given by cq. ©. There is a clear rela- 
tionship. Any level of deep internal energy generation 
increases the central entropy, lifts the degeneracy and ex- 
pands the planetary radius. Figure [1] indicates that the 
radius anomaly, or departure from the zero temperature 
radius, increases with an increasing internal dissipation 
rate that is correlated to the characteristic thermal tidal 
dissipation rate LthT- That is, figure Q] lends strength to 
the hypothesis that thermal tidal torques power the core 
luminosities of hot Jupiters. 

Figure Q] implies that radius excess A, rather than 
planet radius R p is more sensitive in ascertaining the 
actual internal luminosity of a given planet and is there- 
fore, a more useful parameter. The reason for this is 
straightforward: the mass-radius relation of cold spheres 
(e.g., Zapolsky & Salpeter 1969) indicates that - for a 
fixed composition - radius is approximately independent 
of mass. Thus, the gravitational binding energy and ap- 
proximately, the Fermi energy, vary by nearly two orders 
of magnitude. As a result, the energy and power require- 
ments to lift the degeneracy and inflate the radius are 
far more severe for massive objects. The relatively small 
scatter in figure [I] in comparison to similar figures that 
utilize planet radius (e.g., figure 1 of Demory & Seager 
2011; figure 1 of Lauglin et al. 2011; figure 1 of Fortney 
et al. 2011) is partially due to this fact. 

4. DISCUSSION 

Figure Q] is important irrespective of whether or not it 
supports the thermal tide model of Arras & Socrates. It 
indicates that the internal luminosity of the hot Jupiters 
varies by ~ 4 — 5 orders of magnitude. Such a con- 
clusion can be indirectly inferred from e.g., recent work 
by Spiegel & Burrows (2013), independent of the actual 
source of thermal power at great depth. 

With the exception of the thermal tide scenario advo- 
cated here, nearly every plausible model of steady-state 
hot Jupiter inflation invokes or presumes some fixed frac- 
tion of absorbed stellar energy to be redistributed to 
great depth (e.g., Guillot & Showman 2002; Batygin & 
Stevenson 2010; Youdin & Mitchell 2012; Laughlin et al. 
2012). The stellar flux varies like 1/D 2 and therefore, 
so does - for fixed planet radius and stellar type - the 
available amount of energy absorbed. Largely due to the 
fact that hot Jupiters - by definition - possess orbital 
separations 0.01AU < D < 0.1AU, the spread in the 
available energy in starlight varies by approximately two 
orders of magnitude. Since the observed spread in radii 
requires a corresponding spread in steady state power 
of ~ 4 — 5 orders of magnitude, theories that utilize a 



constant fraction of the absorbed stellar light are then, 
effectively, falsified. 

In general, the amplitude of a tidal response is ex 1/D 3 
and therefore, the energy in the tidal interaction, which 
is proportional to the square of the tidal amplitude, is 
oc 1/D 6 . For the case of thermal tides, the tidal response 
is forced by the stellar radiation flux and is therefore 
ex 1/D 2 . The tidal interaction, which leads to orbital and 
spin evolution, is given by the product of the thermal tide 
quadrupole and the tidal potential and is thus oc 1/D 5 . 
Consequently, the expected spread in thermal tidal power 
for the hot Jupiters is roughly 5 orders of magnitude, 
consistent with what is required to inflate their radii. 

The hot Jupiters sit deep within a steep tidal potential 
(ex 1/D 3 ) that leads to even more rapidly varying tidal 
interaction potentials (oc 1/D 6 for gravitational tides and 
oc 1/D 5 thermal tides). Relatively modest displacements 
in orbital separation leads to relatively large changes in 
the internal tidal luminosity. In a way, the radius excess 
phenomena displayed by the hot Jupiters can be thought 
of as a manifestation of, otherwise degenerate, gaseous 
self-gravitating bodies probing deep and rapidly varying 
interaction potentials. 2 

5. SUMMARY 

An observed relationship between radius excess A and 
a characteristic thermal tidal power L thT for the hot 
Jupiters is given in figure [1] It indicates that, in order 
to produce the observed trend in radius excess A, the 
internal luminosity must vary by ~ 4 — 5 orders of mag- 
nitude for the population of inflated hot Jupiters. This 
empirical requirement lends support to any tidal scenario 
that inflates the hot Jupiters steadily over the age of the 
Galaxy and therefore, lends support to the steady-state 
thermal tide scenario of Arras & Socrates (2009a; 2010). 

Furthermore, the empirical requirement that a large 
spread of inferred internal luminosity is required to ex- 
plain the trend in radius excess A from figure [1] falsi- 
fies theories of hot Jupiter inflation that invoke, presume 
or calculate that a fixed fraction of the absorbed stellar 
power is transferred to great depths. 

Both equilibrium and dynamical thermal tides approx- 
imately yield comparable values for equilibrium thermal 
tidal power and planetary spin rate. The role of thermal 
tides can be viewed as allowing for the steady dissipa- 
tion of tidal energy that is continuously replenished by 
the stellar radiation field as it performs work by mov- 
ing matter across the tidal potential (Arras & Socrates 
2009a; 2010). 
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